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Abstract. Optically-thick envelopes may form following 
the tidal disruption of a star by a massive black hole. Such 
envelopes would reprocess hard radiation from accretion 
close to the black hole into the UV and optical bands 
producing AGN-luminosity flares with duration ~ 1 year. 
We show that due to relativistic effects, the envelopes are 
convective. If convection is efficient, then the structure of 
the envelopes is similar to that described in previous work; 
however, the photospheric radius is shown to be very sensi- 
tive to the luminosity at the envelope base, suggesting that 
either the envelope collapses or the envelope expands to a 
maximum radius at which point a wind may set in. For an 
envelope without winds, we find a maximum photospheric 
radius of ~ 10 16 cm (i.e. minimum effective temperature 
~ 6,000 K). The evolution of the envelopes is described 
based on simple energy arguments. 

Key words: Galaxies: active - Galaxies: nuclei - quasars: 
general - Ultraviolet: galaxies 



1. Introduction 

One likely outcome of the tidal disruption of a star by a 
massive black hole is a bright flare with duration of a few 
months to years (e.g. Rees 1988, Ulmer 1998a; hereafter 
U98). The flares are generally thought to be quite hot. For 
example, the temperature associated with an Eddington 
luminosity emitted from a spherical photosphere at the 
tidal radius is 
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where M§ is the black hole mass in units of a million so- 
lar masses and M* and i?* are the mass and radius of 
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the star. Consequently, the discussions of the spectra and 
observability of flares have focused on hard emission (e.g. 
Sembay & West, 1993) or the extreme, ~ 7.5 magnitude, 
bolometric corrections to optical (e.g. U98). Loeb & Ulmer 
(1997, hereafter LU97) make the interesting suggestion 
that if part of the tidal debris forms an extended enve- 
lope around the black hole, then the light may be largely 
reprocessed down to optical with effective temperatures of 
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In the model, an accretion disk forms close to the black 
hole, and the energy released from the disk is reprocessed 
by the extended envelope which connects in an unknown 
way to the disk. 

In Sect. |^, we show that without knowledge of the base 
luminosity which is provided by disk accretion onto the 
black hole, it is impossible to specify uniquely the result- 
ing envelope structure. In Sect. [|, we discuss the small, 
but important, redshift effects from the Schwarzschild ge- 
ometry and show that the envelopes are convective. The 
system's time scales are investigated in Sect. and using 
simple energy arguments, we make predictions regarding 
the evolution of the envelope. A discussion of these results 
is presented in Sect. 

2. Lack of Closure Relations 

The envelope density profile derived in LU97 (Eq. 3) re- 
sults in a logarithmicly diverging mass, so an artificial 
cutoff of the atmosphere was taken in LU97 (Eq. 8). A re- 
lated problem is that the pressure close to the photosphere 
is not handled in a self-consistent manner. In this section, 
we address these problems and discuss their implications. 

We operate in a Newtonian potential, assume that the 
envelope is radiative, and link the envelope model to the 
atmosphere using the Eddington approximation in which 
the surface temperature, To = (l/2) 1 / 4 T e ff, where T e g 
is the effective temperature. At the photospheric radius, 
''phot, where the optical depth is 2/3, T — T e g. The Ed- 
dington approximation requires that the atmosphere be 
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thin so that the surface is not much larger than r p hot ■ The 
hydrostatic equation is: 
dP _ GMp 
dr r 2 



(4) 



where M, is the mass of the black hole alone. We neglect 
the envelope mass since M 6I1v /Mbh ~ 10~ 6 . The equation 
of radiative transfer is 
dP rad nLp 

i = ~~A 91 I 5 ) 

dr 4ircr z 
where the luminosity, L, and absorption coefficient, k, are 
assumed to be constant throughout the envelope. Dividing 
the previous two equations yields the familiar result: 
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dP dP AncGM 
where Le is the Eddington luminosity. The right side of 
equation |6| is constant, because the mass is dominated 
by the mass of the black hole, the luminosity is wholly 
supplied by the accretion disk at the base of the envelope, 
and k is dominated by electron scattering. Integrating, we 
find 

P ga8 =(P-Po)(l-£), (7) 

where the surface pressure, Po, is, in the Eddington ap- 
proximation, equal to P ra d,o = L/QircrZ^. This surface 
term was omitted in LU97. 

Because pressure goes as the fourth power of tempera- 
ture, throughout most of the envelope the surface pressure 
term in Eq. (7] is unimportant, and 
L 



= PP> 
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where (3 is a constant. It then follows from the equation 
of state of the gas (P gas — pk^T/ p,m^) an( i the radiation 
(Pad = aP 4 /3) that everywhere in the envelope except at 
low optical depth, 
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This relation was assumed in LU97. 

It is now possible to solve the hydrostatic equation 
(Eq. ||) for temperature or, equivalently, density by using 
the approximate relation between P and p (Eq. ^|): 
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where tq is a constant of integration which was neglected 
in LU97. The constant is determined by the requirement 
that at the photosphcric radius, rphot, T = T c g. 

Following LU97, we find a relationship between the 
envelope mass and (3. The envelope mass is found by in- 
tegrating Eq. [n] from the base radius, 

R T = P (Mbh/M ) 1/3 ~ 10 13 cm (12) 



to r p hot (The exact location of base radius is an unknown; 
it is the poorly understood interface between the disk and 



envelope. The base radius, as argued in LU97, should be 
close to the tidal radius, because that is where the dynam- 
ical effects of angular momentum are expected to become 
important.) 

M mv = I p4irr 2 dr (13) 
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An envelope which contains no more than a solar mass 
must be extremely close to the Eddington limit with j3 <^ 
10~ 4 . Eq. [l5| shows that the photospheric radius is expo- 
nentially sensitive to /3, suggesting that very small changes 
in base luminosity, which is controlled by a poorly un- 
derstood accretion in the disk, may create large changes 
in the photospheric radius. For example, in this Newto- 
nian approximation, the steady state radius would have 
to change from 10 14 to 10 15 cm when (1 — L/Lg) changes 
from ~ 8 x 10~ 5 to ~ 5 x 10~ 5 . This result illustrates 
the extreme fine tuning (to three parts in 10 5 ) which was 
required in LU97 to reach the steady state solution. As 
we describe in Sect. ^, the fine tuning required in the 
Schwarzschild case is less, but the level is ~ 2% above 
the local critical luminosity. 

We now estimate the maximum radius at which the 
Eddington approximation is valid. The approximation re- 
quires that the atmosphere be thin, or equivalently, that 
the gas pressure scale height be much less than the radius 
at the photosphere. This limit is also a physical dividing 
line, because for further extended atmospheres, one ex- 
pects high mass loss as the atmosphere is nearly isother- 
mal at r <C 1, so the escape velocity falls faster than the 
thermal velocity. 

The gas pressure scale height at the photosphere can 
be evaluated using Eqs. ^ and || as: 
1 dP 

- as (16) 
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At 

?"phot> the gas pressure scale height must be much 
less than the photospheric radius: 
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where we have written (3 in terms of the envelope mass 
(Eq. 15). Note that result is independent of the black hole 
mass. The effective temperature therefore has the same 
scaling as in LU97: 
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These limiting results are surprisingly close to those of 
LU97 who found r pho t ~2x 10 15 cm and T cS ~ 15, 000 K. 
The difference, of course, is that our results show that 
smaller photospheric radii could exist. The exact photo- 
spheric radius cannot be determined without knowledge 
of the inner luminosity source. 

An additional way to see the wind constraint on the 
radius is to consider r^ 1 in Eq. which enforces the 
condition that T = T c g at the photosphere (which radi- 
ates at the Eddington limit). The constant is zero when 
fphot ~ 10 15 , and becomes negative for larger values of 
fphot- Negative values of r^ 1 give (see Eq. |l(]) regions at 
the top of the envelope which approach an isothermal state 
and which would likely drive winds. 

Even before the radius expands so far that the pressure 
scale height becomes comparable with the radius, winds 
could begin to play an important role-either in altering 
the envelope structure or in removing much the envelope 
mass. A estimate of the importance of the winds may be 
found by connecting, at the photospheric radius, isother- 
mal wind solutions to the envelope solutions. Because of 
the extremely low densities (~ 10~ 15 g/cm 3 ), the enve- 
lope is ionized and the cross-section is largely provided 
by electron-scattering, so an isothermal wind powered by 
the continuum cross-section (rather than line transitions), 
may be most appropriate. In this case, the following equa- 
tion describes the outflow (e.g. Kudritzki, 1988): 

- 2 N dv 2v 2 GM b „„ L 
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v dr r r 2 v Le 

where v s is the sound speed (in our case, isothermal so 
that v s — (kBTcg/iirrih) 1 / 2 and the sonic point, r s = 
GM(1 — L/L E )/2v 2 . Dimensionless radius, ij — r/r s , and 
dimensionless velocity, u = v/v s , allow Eq. |T]to be neatly 
integrated to obtain: 

■ 2 2 

21nry + 1.5. (22) 
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Connecting such outflows to the Eddington envelopes al- 
lows for determination of mass loss as a function of en- 
velope photospheric radius. Generally, the time for such 
a wind to significantly reduce the envelope mass is at 
least 1000 years — many times longer than the other rel- 
evant time scales for the problem (results are shown in 
Fig. 1). For photospheric radii larger than ~ 3 x 10 15 cm, 
the sonic point occurs below the photosphere, showing, in 
agreement with the limiting radius found above, that the 
hydrostatic solution (Eqs. |l0| and [n]) is no longer valid. 



3. Effects in the Schwarzschild Metric 

The Schwarzschild geometry has a significant impact on 
envelopes which radiate near the Eddington limit. At first 
appearance, the Schwarzschild geometry seems to be an 
unimportant correction to the problem, because the base 
of the photosphere is at ~ 25i?Si and redshift effects are 



of order 2%, (1 + z) ~ (1 + Rs/2r). However, because the 
envelopes are very close to the Eddington limit, as is re- 
quired to produce an extended photosphere, the envelope 
structure is very sensitive to /3 = 1 — L/ L E . In particular 
for envelopes of ~ 1M©, /? is ^ 10~ 4 , so a 2% reduction in 
luminosity will significantly alter the envelope structure. 

We quote below a number of results which were ob- 
tained in the study of x-ray bursts and envelopes around 
neutron stars (Paczyhski & Anderson 1986; Paczyhski & 
Proszynski 1986). For our purposes, the most important 
feature of the relativistic stellar structure equations is that 
the local critical luminosity does not scale with radius in 
the same manner as the local luminosity with the conse- 
quence that the structure is convective rather than radia- 
tive. Specifically, the local luminosity and critical luminos- 
ity determined by the local gravitational forces (equivalent 
to the Eddington luminosity at large r) are 

L = L x ! .1 - — ! (23) 
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where is the luminosity as measured far from the black 
hole. Because of the different scalings, a near critical lumi- 
nosity at the surface requires a super-critical luminosity 
at the base. As a consequence the star is convective. Alter- 
natively, the convective nature can be seen by comparing 
the radiative and adiabatic gradients: 
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Even for a near critical surface luminosity, V ra d quickly 
becomes larger than V a d as the coordinate, r, decreases. 

As a consequence, the envelope must be convective. 
Whether or not the convection is efficient is difficult to 
determine because in radiation dominated regimes, con- 
vection is not fully understood, although progress is being 
made (e.g., Arons 1992). If convection is efficient, then by 
definition, V a( j is nearly a constant. The equation of the 
temperature gradient yields: 



dP ra ,d 
dP 



1 
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j/3 32 /3- 

so P TEl d/P is very closely a constant (to 



32 M 



(29) 



1 part in 10 -9 ), 

and we recover the polytropic equation of state P = Ap 4 ^ 3 , 
which was used in the previous section. Assuming that 
convection is efficient, we can determine the relationship 
between the photospheric radius, the envelope mass, and 
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the surface luminosity in the same manner as in Sect. 0. 
Ignoring the relativistic terms in the hydrostatic stellar 
structure equations, which is a good approximation (to 
~ 2%) because the terms enter multiplicatively rather 
than as differences, we recover eqn. [l5[ In contrast to the 
Newtonian case, in the Schwarzschild case, the luminosity 
at the base must be locally super-critical in order to sup- 
port an extended envelope with a near-Eddington surface 
luminosity. 

If the luminosity is ~ 2% super-critical, and the en- 
velope is able to expand, then photospheric radius is still 
sensitive to the base luminosity, but not exponentially sen- 
sitive as is the case for a Newtonian atmosphere. The lu- 
minosity at the photosphere is nearly equal (to better than 
one part in 10 4 ) to the critical luminosity (Eqs. 24), so 

/ d \-i/2 
iphot « L CI = L E 1 2_ . (30) 

\ ^phot / 

Using equation |25| to write the luminosity at the base, Lb 
as a function of photospheric radius, yields the result: 

Rs 



VotA 1000 years) 



fphot 



1 



(31) 



(W^b) 2 ' 

In steady state the photospheric radius would increase 
from 10 14 to 10 15 cm with a 0.1% change in base lumi- 
nosity. 

There may be additional effects which compete with 
the relativistic effects. For neutron stars, the temperatures 
at the base of the envelope are so high that relativistic cor- 
rections to the Thompson cross section become important 
(e.g. Paczyhski & Anderson 1986), but the temperatures 
around the tidal disruption created envelopes never reach 
such high temperatures. In our case, slight rotation may 
serve to reduce the critical luminosity at small radii, in 
the plane of rotation. Along the rotation axis, the enve- 
lope would likely be convective for the reasons described 
above. 



4. Evolution in Time 

Both the Newtonian and relativistic envelopes require fine 
tuning of the luminosity in order to produce static ex- 
tended envelopes, because r p hot depends sensitively on the 
luminosity. In the Schwarzschild case, unless the base of 
the envelope knows to shine locally at ~ Rs/2r\ i ~ 2% 
above the local critical luminosity, and the luminosity is 
restricted to a narrow regime, it appears difficult to main- 
tain the type of envelopes discussed in Sect. 0. The de- 
pendence on luminosity is strong; a fraction of a percent 
change in base luminosity could increase the steady-state 
radius from 10 14 to 10 15 cm. Since the energy source (the 
accreting torus) is probably not extremely well-coupled to 
the envelope, it is unlikely that the base luminosity is so 
finely tuned. Therefore, we consider how the photosphere 
will change with time for different base luminosities. 




(cm) ] 



Fig. 1. Time scales for envelopes as a function of pho- 
tospheric radius: dynamical time, thermal time, and time 
for wind driven mass loss. The thermal time scale is calcu- 
lated assuming a energy input of 2Le and energy output 
of Le and rb = 10 13 cm. The envelopes can stay in hydro- 
static equilibrium only as long as idyn < hherm- The mass 
loss to winds produced by an isothermal outflow are small 
enough that the do not appear to affect the system. 



We identify three cases. If the base luminosity is sub- 
critical, the envelope is not much bigger than the accreting 
torus. When the luminosity is super-critical, the envelope 
must expand. If the supcr-criticality is less than ~ 2%, 
the expansion is modest and the envelope may remain in 
hydrostatic and thermal equilibrium. If the base luminos- 
ity is super-critical by a significant factor, e.g. 2, then the 
envelope expands hydrostatically on a thermal time scale, 
until it expands so much that the dynamical time scale 
becomes shorter than the thermal time scale. After that, 
the outer parts of the envelope are no longer in hydro- 
static equilibrium, but this does not imply that they are 
instantly lost. The envelope may continue to expand, per- 
haps even with a structure similar to the hydrostatic so- 
lution, until the scale height of atmosphere becomes com- 
parable to the radius (as discussed in Sect. ||) at which 
point a wind will set in. 

We now examine the expansion of the envelope. The 
relevant time scales are the dynamical time scale, the ther- 
mal time scale, and the radiation time scale, over which 
time all matter would be accreted. The first two time 
scales are shown in Fig. 1 as a function of photospheric 
radius. 
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The dynamical time scale grows with photospheric ra- 



dius: 



^dyn 



GM BH (1-L/L E ) 
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years. (33) 

The thermal time scale, defined as the time scale to 
unbind the envelope assuming the energy in the envelope 
increases at a rate of L + = L in — L out w Le (e.g. if the 
base luminosity were 2Le and the surface radiated at Le)- 
The total energy of the hydrostatic envelope is E tot = 
-Egrav + ^internal, where the internal energy is 

^internal ~ / AlTT 2 P rad dr (34) 

3 GM B nM cnv 
W 4r b [ln(r phot /r b )-1.8]- W 

Note that the energy is strongly dependent on the loca- 
tion of the inner boundary, r D . The gravitational energy, 
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with a complex dependence on photospheric radius and 
black hole mass. The thermal time scale depends both 
on the net energy injection rate and on the base radius of 
the envelope. As both parameters depend on the unknown 
physics of the torus, and as such, their exact values are 
uncertain. We believe the base radius to be near the tidal 
radius ~ 10 13 cm, where rotational support likely becomes 
important. Similarly, thick tori may produce luminosities 
up to a few times the Eddington limit, but the exact value 
cannot be predicted. These parameters enter into the ther- 
mal time scale multiplicatively, and together may lengthen 
or shorten the thermal time scale by a factor of 2 or 3. 

Lastly, the radiation time is the time to radiate, at the 
Eddington limit, all of the accretion energy of the tidal 
debris: 

i rad » 20AV 1 years, (37) 

for an accretion efficiency of 10% and an envelope mass 
of 0.5M Q . This time scale does not depend on the pho- 
tospheric radius, and is generally longer than the other 
relevant time scales of interest (for a 10 6 M Q black hole). 

As discussed above, if the luminosity at the base is 
strongly super-critical, then the envelope will expand and 
the photosphere will move to larger radii (see Fig. 2). For 
an envelope with its photosphere below ~ 10 14 cm, the 
dynamical time is shorter than the thermal time, so the 
envelope should expand in hydrostatic equilibrium. When 
the photosphere grows larger than ~ 10 14 cm, the enve- 
lope can no longer be in hydrostatic equilibrium. However, 
the material at small radii has a much shorter dynamical 
time, so may be able to adjust to a near-equilibrium state 
quickly. The envelope could continue to expand, perhaps 
with the interior in a near equilibrium state and the outer 
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Fig. 2. Radius and temperature for sequences of expand- 
ing envelopes labeled with three different energy input 
rates (each envelope radiates at the Eddington limit) and 
a base radius of 10 13 cm. The solid lines are the phase in 
which the envelopes are in hydrostatic equilibrium. The 
dashed lines signify that the envelope is not in hydro- 
static equilibrium, but for the purposes of this illustra- 
tion, we assume that the models do stay in hydrostatic 
equilibrium and calculate the time between models at 
At = AEy(L in _L|rmout), where AE is the energy differ- 
ence between models. These calculations are described in 
more detail in Sect. H. 



parts further from equilibrium. If the photospheric radius 
reaches ~ 10 15 cm, the thermal time scale is so short and 
the pressure scale height becomes so large compared to 
the radius, that it seems most likely either a strong wind 
would form and carry the excess energy away, or the entire 
envelope would be unbound. 

Fig. 2 illustrates possible evolutionary sequences for 
different values of the base luminosity. For Fig. 2, we make 
the assumption that the envelope is in a hydrostatic state 
(described by eqn. O), even after the thermal time scale 
falls below the dynamical time scale. The time to move 
between envelopes At — AE/L + , where AE is the en- 
ergy difference between sequential envelopes (A similar 
method was applied by Ulmer (1998b) to investigate the 
evolution of thick accretion disks.) Evolutionary sequences 
show the evolution of both photospheric radius and tem- 
perature with time. The curves become nearly vertical as 
the envelopes evolve to large radii, because the binding 
energy becomes so small that they can be unbound by 
very little energy input. 



6 



Ulmer, Paczyriski, & Goodman: Tidal Disruption Eddington Envelopes 



It seems likely that the envelope will expand on a time 
scale of months to years, depending on the base radius and 
the base luminosity, to ~ 10 14 cm where T c s ~ 50, 000 K. 
Subsequently, the envelope will expand, but will not be 
able to maintain hydrostatic equilibrium throughout. 



5. Discussion 



Eddington envelopes require fine tuning of the luminos- 
ity in order to produce extended envelopes. It appears 
unlikely that the luminosity is held within such a tight re- 
gion. More reasonably, either the average luminosity will 
be sub-critical at the base, and the envelope will collapse, 
or the luminosity will be super-critical at the base, and 
the envelope will grow (and thereby radiate at lower tem- 
peratures between ^80,000 K and ^20,000 K as shown 
in Fig. 2) and eventually reach the most extended hydro- 
static model or the wind solutions, as described in Sect. g. 
This scenario may be valid only if convection is efficient. 
The stellar structure equations are not as tractable if con- 
vection is inefficient, and we have not attempted to solve 
them here. 

We reiterate that the remarkable feature of extended 
envelopes is that they could reprocess much of the hard 
radiation produced by a tidal disruption event into the 
optical band. Bolometric corrections would drop from ~ 
7.5 magnitudes to 1-2. The color temperatures of the 
objects may be slightly higher than the effective tem- 
peratures due to suppression of absorption processes at 
such low densities, but such changes are expected to be 
about a factor of two (see the work on neutron stars by 
Babul & Paczyhski (1987)). Even scattering-dominated 
envelopes retain their important property of producing 
optically bright flares. 
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